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A  FORTRAN  computer  code  (SKINTEMP)  har  been  developed  to  calculate 
transient  missile/aircraft  aerodynamic  heating  parameters  utilizing  basic  flight 
parameters  such  as  altitude,  Mach  number,  and  angle  of  attack.  The  insulated  skin 
temperature  of  a  vehicle  surface  on  either  the  fuselage  (axisymmetric  body)  or 
wing  (two-dimensional  body)  is  computed  from  a  bas*''  ’leat  balance  relationship 
throughout  the  entire  spectrum  (subsonic,  transonic,  jv  '  sonic,  hypersonic)  of 
flight.  This  calculation  method  employs  a  simple  finite  difference  procedure  which 
considers  radiation,  forced  convection,  and  non-reactive  chemistry.  Surface 
pressure  estimates  are  based  on  a  modified  Newtonian  flow  model.  Eckert’s 
reference  temperature  method  is  used  as  the  forced  convection  heat  transfer  model 
SKINTEMP  predictions  are  compared  with  a  limited  number  of  test  cases. 
SKINTEMP  was  developed  as  a  tool  to  enhance  the  conceptual  design  process  of 
high  speed  missiles  and  aircraft.  Recommendations  are  made  for  possible  future 
development  of  SKINTEMP  to  further  support  the  design  process. 
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I .  INTRODUCTION 


A.  GENERAL 

A  FORTRAN  computer  program  has  been  created  to  determine 
transient  missile  and  aircraft  mission  aerodynamic  heating 
parameters  such  as  heat  transfer  coefficients,  insulated  skin 
temperatures  and  heat  flow  rates  at  some  point  on  the  surface 
a  wing  or  a  body  of  revolution.  Required  data  supplied  from 
the  user  include  material  properties  and  trajectory  data  such 
as  Mach  nixmber,  angle  of  attack,  altitude,  and  time  of  flight. 
This  program  calculates  aerodynamic  heating  parameters  for 
sxibsonic,  transonic,  supersonic  and  hypersonic  (no  reacting 
chemistry)  flight  regimes. 

B.  BASIC  THEORY 

The  governing  equation  which  serves  as  the  basis  for  the 
computer  program  incorporates  both  radiative  and  forced 
convective  heat  transfer  processes.  The  thermal  energy 
balance  equation  for  this  insulated  skin  heat  transfer  case  is 
given  by 

=  hA(r,^  -  T^)  -  06A(r/  -  T/)  .  (1) 

The  equation  written  in  this  form  assumes  that  the  adiabatic 
wall  temperature,  T,^  is  greater  than  the  insulated  skin 
temperature,  T,,  which  is  greater  than  the  radiation  reference 


1 


temperature,  T,.  The  heat  transfer  coefficient  for  convection 
is  obtained  from  the  Nusselt  nxnnber  which  is  a  function  of  the 
Reynolds  and  Prandtl  numbers.  These  dimensionless  numbers  are 
defined  as 

(a)  (i)  (c)  (2) 

Once  the  complex  heat  transfer  coefficient  calculation  is 
performed  and  the  other  significant  heat  transfer  parameters 
determined,  a  finite  difference  numerical  method  can  be 
applied,  yielding  a  skin  temperature  profile  for  the  specified 
point  of  interest. 

C.  SCOPE 

Chapter  II  details  the  theory  used  to  develop  the  SKINTEMP 
program. 

Chapter  III  describes  the  actual  coitputer  program,  its 
required  entry  data,  expected  output  data  and  general  program 
limitations . 

Chapter  IV  presents  test  case  results  with  the  appropriate 
comparisons,  which  include  intermediate  heat  transfer 
calculations.  One  test  case  is  obtained  from  a  heat  transfer 
textbook  and  compared  to  corresponding  SKINTEMP  results.  A 
less  descriptive  NASA  test  case  is  also  compared  with  SKINTEMP 
results . 

Chapter  V  provides  recommendations  for  future  code 
development . 
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II.  TRANSIENT  HEAT  FLOW  THEORY 


The  fundamental  concepts  used  to  build  the  theoretical 
model  for  SKINTEMP  involve  two  basic  modes  of  heat  transfer, 
convection  and  radiation.  Forced  convection  is  assumed  because 
the  vehicle  is  propelled  through  the  air  by  the  release  of 
chemical  energy.  During  supersonic  and  hypersonic  flight, 
the  local  stagnation  pressure  at  the  edge  of  the  boundary 
layer  is  assumed  to  be  defined  by  the  stagnation  pressure 
behind  a  normal  shock  wave.  Figure  1  illustrates  this 
concept . 
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Boundary  layer  flow  is  assumed  to  be  turbulent  due  to  the 
magnitude  of  the  Reynolds  number  which  is  generally  greater 
than  500,000.  The  gas  dynamic  relations  are  based  on  inviscid 
flow,  no  reacting  chemistry  assumptions  [Ref.  1] . 

The  primary  objective  of  the  computer  code  is  to  rapidly 
compute  the  transient  insulated  skin  temperatures  along  the 
trajectory  of  the  aircraft  at  a  given  point  somewhere  on  the 
surface.  To  accomplish  this,  the  thermal  energy  balance 
equation 

=  hA{T^^-  T^)  -  aeA{T/  -  T/)  (l) 

is  solved  for  T,  using  a  finite  difference  method  [Ref.  2] . 
However,  before  this  can  be  accomplished,  the  adiabatic  wall 
temperature  and  the  heat  transfer  coefficient  must  be 
determined.  The  required  data  for  the  skin  temperature 
calculations  include  trajectory  parameters  (time,  Mach 
number,  angle  of  attack,  altitude) ,  material  propercies 
(density,  material  specific  heat,  emissivity) ,  the  location 
and  skin  thickness  at  the  point  of  interest,  and  a  radiation 
reference  temperature.  This  radiation  temperature  is  the 
temperature  of  a  distant  surface  seen  by  the  insulated  skin 
element,  which  for  this  program,  is  either  space  or  the  earth. 
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A.  THE  ADIABATIC  WALL  TEMPERATURE 

The  first  parameter  needed  to  determine  skin  temperature 
is  the  adiabatic  wall  temperature.  The  adiadjatic  wall 
temperature  is  given  by 


I'm  =  n  [l  *  ■ 

The  local  Mach  nxjmber  needed  for  the  solution. 


(3) 


Mr  = 


-1 


(4) 


Y-1 


is  dependent  on  the  local  pressure  and  the  local  stagnation 
pressure.  The  speed  of  the  aircraft,  whether  subsonic  or 
supersonic  will  determine  the  local  stagnation  pressure.  If 
the  trajectory  profile  calls  for  subsonic  flight,  then  the 
local  stagnation  pressure  will  be 


PtL  =  Pl 


1  +  ^^M^ 


T/T-l 


(5) 


If  the  Mach  number  is  in  the  supersonic  range,  the  local 
stagnation  pressure  is  then 


where 


tL 


=  P. 


1  + 

•  2 


T/T-l 


111 


1/1-T 


Y+lJ 


(6) 


i’t-  =  P- 


1  + 


y/t-1 


(7) 
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Freestream  pressure,  temperature,  and  density  are  calculated 
using  standard  atmosphere  property  value  approximations  found 
in  Reference  3 . 


Local  surface  pressure  is  determined  using  modified 

Newtonian  theory,  [Ref.  4]  where 

Pl  =  C^cos^Qg  *  P,  (3) 

_  ^CL  ~  -P- 


^pmax 


g  =  l/2p.vf 


(9) 

(10) 


Modified  Newtonian  flow  theory  has  been  shown  to  be  applicable 
for  the  prediction  of  local  surface  static  pressures  over  all 
surfaces  experiencing  non-separated  flow.  [Ref.  5,6]  The 
accuracy  of  Modified  Newtonian  flow  theory  is  quite  reliable 
at  hypersonic  speeds.  Its  use  at  subsonic  and  low  supersonic 
speeds  is  justified  by  the  continuity  and  simplicity  it 
provides  in  the  lower  Mach  number  regions  where  aerodynamic 
heating  rates  are  so  low  as  to  be  negligible. 

The  remaining  parameter  needed  to  determine  the  adiabatic 
wall  temperature  is  the  recovery  factor  (turbulent  flow  is 
assumed) ,  given  by 

r  =  (11) 


where  the  Prandtl  number  is 


(12) 

The  specific  heat  of  air  used  in  the  Prandtl  number  is  given 
by  [Ref.  2] 

Cp  =  .219  +  .  342^-4 -  .  293^-8 T^.  (13) 
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Thermal  conductivity  [Ref.  7]  and  dynamic  viscosity  [Ref.  2] 
are  also  temperature  dependent  and  are  given  rr iipectively  by 


k  =  .0141 


7  07 

r^+225 

492, 

13/2 


H  = 


2.21  E-Q 
r^+198.7 


B.  TEE  HEAT  TRANSFER  COEFFICIENT 


(14) 

(15) 


In  order  to  calculate  the  heat  transfer  coef f icient , one 
uses  the  concept  of  the  reference  temperature.  The  reference 
temperature  is  a  temperature  whose  magnitude  is  typically 
between  that  of  the  wall  temperature  and  freestream 
temperature  and  in  low  speed  flow,  the  reference  temperature 
does  become  the  average  of  the  two  [Ref.  2] .  This  reference 
temperature  more  accurately  models  the  temperature 
distribution  that  occurs  within  the  boundary  layer,  which  has 
a  significant  impact  on  the  skin  temperature  [Ref.  8]  .  The 
transport  properties  should  be  evaluated  at  the  reference 
temperature  [Ref.  9]  which  is  defined  by 


T*  =  +  .5(Tg  -  rj  +  .22{Tj^„-  T^)  .  (16) 


Values  for  the  thermal  conductivity,  dynamic  viscosity,  and 
specific  heat  used  in  the  Nusselt  number  are  based  on  this 
reference  temperature. 


The  Nusselt  number  can  be  defined  as 


Based  on  experimental  data,  the  Nusselt  number  can  also  be 
written  as  a  function  of  Reynolds  and  Prandtl  number,  [Ref.  2] 

=  .0292(i^jJg)  (18) 

where  the  boundary  layer  is  assumed  to  be  turbulent  and  the 
Reynolds  number  is  defined  as 


(19) 


With  these  significant  variables  defined,  the  thermal 
energy  balance  equation  can  be  solved  using  a  simplified 
finite  difference  method.  This  method  calculates  the 
temperature  value  at  a  future  time  based  on  a  calculated  time 
increment  and  the  current  temperature  value. 


C.  THE  SKIN  TEMPERATURE 

Solution  of  the  thermal  energy  balance  equation  for  the 
incremental  skin  temperature  yields 

-  T*)]^  +  (20) 


For  stability  purposes,  the  time  interval  associated  with  this 
finite  difference  method  is  given  by 


At  ^ ^ 


h+h. 


(21) 


where 


=  Aea  [rj  +  2’^] 


(22) 


Each  progressive  skin  temperature  is  dependent  on  the 
previous  temperature,  heat  transfer  coefficient  and  adiabatic 
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wall  temperature.  The  time  interval  choice  is  critical  and 
must  meet  the  stability  requirement  of  the  finite  difference 
method.  This  stability  requirement,  which  prevents  the 
solution  from  diverging,  is  given  by 


At(h  +  H^) 


1. 


(23) 


To  meet  this  requirement,  the  time  interval  was  chosen  as 


At 


.25 


mCp 

h  +  H/ 


(24) 


D.  THE  HEAT  FLOW  RATE 


The  heat  flow  rate  is  dependent  on  the  location  of  the 
stagnation  heating  point.  In  the  stagnation  region  for  the 
axisymmetrical  and  two  dimensional  platforms,  the  maximum 
stagnation  point  is  70  degrees  from  the  centerline  [Ref.  10] 
and  the  stagnation  heating  rate  is  given  by  [Ref.  11] 


*stag 


_  17600 

_e. 

u 

3.15 

[  "si  -  "vs 

y/Rad  > 

Po 

[“cj 

"si  ~  ^W545, 

(25) 


Away  from  the  stagnation  point  but  prior  to  regions  that  can 
be  treated  as  flat  plates,  the  heat  flow  rate  [Ref.  10]  is 
assumed  to  be 


Qheat  =  g^tagCOS(0)3/^ 


(26) 


otherwise,  in  areas  where  flat  plate  assumptions  hold,  the 
heat  flow  rate  is 

-  oeAiT/  -  r/)  .  (27) 
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III.  THE  SKZMTBKP  COMPUTER  PROGRAM 


A.  THE  COMPUTER  CODE  STRUCTURE 

The  flow  logic  used  for  the  SKINTEMP  program  is 
illustrated  in  Figure  2.  Th''  program  is  divided  into  four 
major  sections,  data  entry,  atmospheric  property  generation, 
skin  temperature  calculations  and  data  output. 

The  first  module  contains  the  data  entry  loop.  The  data 
can  be  entered  via  two  methods,  keyboard  entry  or  file  entry. 
If  the  user  chooses  to  utilize  keyboard  entry,  the  program 
will  prompt  the  user  for  the  specific  data.  Each  entry  should 
be  followed  by  a  keyboard  hard  return.  The  file  entry  method 
requires  the  user  to  create  an  external  file  using  any 
standard  editor  containing  all  the  required  information.  For 
file  entry,  the  correct  syntax  to  enter  at  the  pron^t  is 
SKINTEMP. in.  The  .in  refers  to  the  entry  file;  any  input  file 
used  must  have  the  .in  included  in  the  file  name.  This 
command  allows  data  to  be  read  automatically  by  the  input 
subroutine. 

The  important  note  to  remember  is  that  for  file  entry, 
each  piece  of  fixed  point  data  must  be  placed  on  a  separate 
line  cind  the  entry  data  must  correlate  with  the  data  requested 
by  the  program.  If  the  program  terminates  with  an  error 
message,  this  would  be  the  first  place  to  check  for  mistakes. 
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The  second  section  of  the  code  generates  standard 
atmospheric  values  that  are  needed  for  the  freestream 
pressure,  temperature,  and  density  calculations.  The 
governing  equations  used  to  model  atmospheric  properties  are 
divided  into  two  regions,  gradient  and  isothermal  [Ref.  3]. 
The  gradient  relationships  are 


_g_  _ 

Po 

p  _  T  ~  ^3/LR) 
"p”  ~  "t” 

■‘o 

T  =  L{h-h^)  +  T^. 


(28) 

(29) 

(30) 


There  are  several  isothermal  regions  in  the  atmosphere. 
While  the  temperatures  remain  constant  in  each  of  these 
isothermal  regions,  the  tenperatures  of  the  gradient  regions 
will  vary.  [Ref.  3]  The  pressure  and  density  relationships 
for  isothermal  regions  are  given  by 

ir  ®  (31) 

<32) 

Po 

with  these  relationships,  the  atmosphere  is  modeled  to  260,000 
feet  with  isothermal  layers  from  approximately  36,000  to 
82,000  feet  ( teir^erature  =  392.4  ‘Tl,  density  =  .00071  slug/ft\ 
pressure  =  476  Ib/ft^)  and  154,000  to  174,000  feet 
(temperature  =  509  ®R,  density  =  .00000303  slug/ft^,  pressure 
=  2.6419  Ib/ft^)  . 
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The  third  major  section  of  SKINTEMP  is  the  calculation  of 
the  skin  temperature.  This  section  is  a  large  counted  loop 
because  each  successive  skin  temperature  is  dependent  on  the 
previous  temperature  calculated.  This  loop  is  partitioned 
into  four  sections.  The  first  section  calculates  adiabatic 
wall  temperature  and  the  second  determines  the  heat  transfer 
coefficient.  These  two  calculations  lead  to  the  skin 

temperature  solution,  which  is  in  the  third  section.  The 
fourth  and  final  section  calculates  the  heat  flow  rate.  [Ref. 
Chapter  II]  Appendix  A  contains  the  actual  computer  program. 
The  initial  documentation  of  the  program  contains  a  list  of 
the  variables,  their  definitions,  and  the  English  units  used. 

It  is  important  that  the  user  understand  the  differences 
between  the  time  values  used  throughout  SKINTEMP.  The  user 
should  input  time  intervals  that  best  linearize  the 
trajectory.  In  other  words,  whenever  the  Mach  number, 

altitude  or  angle  of  attack  is  changing,  this  time  and 
associated  flight  parameters  should  be  a  data  entry  point. 
These  time  intervals  can  be  the  print  out  time  intervals  but 
SKINTEMP  also  allows  the  user  to  print  the  results  at 
different  time  intervals  if  desired.  These  input  and  output 
time  intervals  are  independent  of  the  time  interval  chosen  for 
the  heat  balance  calculations.  Figure  3  demonstrates  these 
differences  and  should  be  used  as  a  reference  when  developing 
a  data  entry  strategy  for  SKINTEMP. 
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Figure  2 .  SKINTEMP  Flowchart 


See  Chapter  II  for  a  complete  description  of  the  relationships 
needed  to  obtain  the  desired  results. 
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^  TIME  INTERVAL  COMPARISONS 


INPUT  TIME 

H — H — I  -I  -(  H  -41 - M — H — f— 

STEP  TIME 


OUTPUT  TIME 


Figure  3.  Time  Interval  Comparison 


14 


B.  REQUIRED  DATA  INPUT 

SKINTEMP  allows  the  user  to  enter  data  through 
the  keyboard  or  through  an  external  file.  In  either  case,  the 
following  data  is  required  to  successfully  run  the  program: 

•  The  shape  of  the  body,  which  must  be  either  two 
dimensional  or  axisymmetric 

•  If  a  two  dimensional  wing  shape  is  selected,  the  wing 
sweep  (degrees) 

•  The  launch  surface  skin  temperature  of  the  vehicle  (°F) 


•  For  each  time  interval,  the  Mach  number,  the  angle  of 
attack  (degrees) ,  and  the  altitude  (feet) 

WARNING:  Entering  a  Mach  number  of 
zero  will  cause  the 
program  to  abort.  See 
Section  D.  below. 

•  The  distance  from  the  nose  or  leading  edge  to  the  point  of 
interest  (feet) 

•  The  thickness  of  the  surface  at  the  point  of  interest 
( inches ) 

•  The  nose  radius  (feet) 

•  The  material  density  (pound  per  cubic  foot) 

•  The  material  emissivity  (dimensionless) 

•  The  material  specific  heat  (BTU  per  pound-"F) 

•  The  orientation  of  the  surface.  Does  the  surface  face  the 
earth  or  space  for  the  majority  of  the  flight?  This 
provides  the  radiation  reference  temperature  of  -460 "F  for 
space  and  60  "F  for  the  earth.  (°F) 

•  The  angle  from  the  unit  surface  normal  of  the  skin  element 
to  the  freestream  (degrees) 

•  The  output  time  interval  (seconds) 
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C.  SKIKTENP  OUTPUT 


The  output  of  the  computer  prograun  is  in  tcUaular  form. 
The  user  can  choose  any  time  interval  for  printing  independent 
of  the  data  time  entries.  See  Appendix  B  for  an  sairple 
output.  SKINTEMP  will  also  produce  an  external  UNIX  file 
named  Fort. 3  which  produces  a  table  of  values  for  the  time  and 
skin  temperature  that  can  be  used  in  typical  engineering 
graphics/plotting  programs. 

The  program  is  clearly  documented  and  the  print  subroutine 
is  the  last  section  in  the  program.  The  output  format  is  as 
follows; 

TIME  ADIABATIC  WALL  HEAT  TRANSFER  SKIN  HEAT  FLOW 

TEMPERATURE  COEFFICIENT  TEMPERATURE  RATE 

(sec)  (®Rankine)  (BTU/sec-ft^-®F)  ("Rankine)  (BTU/sec) 

D.  SKINTEMP  LIMITATIGNS 

Th3  purpose  of  this  code  is  to  provide  a  design  tool  for 
conceptual  missile  and  aircraft  design  projects.  At  the 
present  time,  SKINTEMP  can  only  calculate  the  transient  heat 
parameters  for  bodies  that  are  either  two  dimensional  or 
axisymmetric  body  shaped.  See  Chapter  V  for  future 
recommendations  for  code  expansion. 

Another  limitation  and  perhaps  the  most  important  for  user 
information  is  that  at  no  time  can  the  Mach  number  entered  be 
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equal  to  zero.  The  reason  for  this  is  that  in  the  calculation 
for  Cp .  dynamic  pressure  is  the  divisor.  Cp„..  and  dynamic 
pressure  are  given  by 

C’p.max  =  ^PtL  -  -P-)/*?  (9)  wheze  q  =  l/2pV^  = 

A  Mach  number  of  zero  will  cause  Cp^  to  be  undefined.  If 
this  occurs,  the  program  will  terminate  without  results.  If 
a  Mach  number  of  zero  is  the  actual  value  achieved  in  the 
trajectory,  a  small  value  should  be  entered,  i.e.,  .001  Mach. 

The  last  limitation  is  the  array  size  of  the  input  data. 
At  the  present  time,  the  array  size  of  data  points  entered  is 
100.  If  more  than  100  data  points  are  desired  by  the  user, 
the  array  size  should  be  altered  in  the  body  of  SKINTEMP.  The 
array  dimension  statement  is  the  first  to  occur  after  the 
initial  documentation  and  the  mcocimum  limit  is  dictated  by  the 
version  of  FORTRAN  being  used. 
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IV.  TEST  CASE  RESULTS 


A.  TRUITT  TEST  CASE 

The  test  case  used  to  validate  SKINTEMP  is  based  on  an 
example  skin  temperature  problem  from  Robert  W.  Truitt's  text 
Fundamentals  of  Aerodynamic  Heating,  [Ref  12]  .  The  example 
missile  profiled  is  axisymmetric  in  shape  and  40  feet  in 
length.  The  skin  material  used  is  Inconel  steel  and  the  skin 
thickness  of  the  material  at  the  point  of  interest,  20  feet 
from  the  nose,  is  .00504  inches.  The  flight  trajectory  data 
are  given  in  Figure  4 . 

Trajectory  Mach  number 


Trajectory  Altitude 


Figure  4.  Truitt's  trajectory  data 
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Truitt  assumes  turbulent  flow  but  does  not  use  specific 
local  or  reference  state  conditions.  Instead,  he  assvunes 
freestream  conditions  throughout  the  calculations.  This 
approach  is  responsible  for  at  least  some  of  the  noticeable 
deviation  in  the  comparative  results  [Ref.  Section  2]. 

A  Runge  -  Kutta  numerical  integration  method  for  the  skin 
temperature  calculation  is  employed  by  Truitt  rather  than  the 
finite  difference  method  used  in  SKINTEMP.  This  difference  in 
approach  was  one  of  the  primary  reasons  this  test  case  was 
used.  Similar  results  from  two  different  equivalent  numerical 
methods  tend  to  confirm  results. 

1 .  Data  Comparisons 

Throughout  Truitt's  skin  temperature  calculation 
example  several  intermediate  heat  transfer  and  fluid  flow 
parameter  values  are  defined.  These  computed  values  include 
freestream  static  pressure  and  temperature,  Reynolds  nvimber 
and  the  computed  skin  temperature  (based  on  freestream 
properties) .  The  comparisons  of  these  parameters  are  listed 
in  the  following  graphs  and  tables  to  validate  SKINTEMP  and 
illustrate  some  of  the  theoretical  differences  between 
calculation  methods . 
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a.  Frees tream  Pressure  Conparisoa 

Figure  5  graphically  displays  the  difference  in 
freestream  static  pressures  while  Table  1  details  the  exact 
values  computed.  Truitt's  values  for  freestream  static 
atmospheric  density  are  not  as  close  to  the  NACA  standard 
atmosphere  [Ref.  13]  as  the  SKINTEMP  code  calculations  and  the 
slight  divergence  in  values  contributes  to  the  difference  in 
the  predicted  skin  temperatures. 


FREESTREAM  PRESSURE  COMPARISON 


S  10  IS  20  25  30 

Time  (seconds) 


Figure  5.  Freestream  static  pressure  comparison 


TABLE  1.  FREESTREAM  PRESSURE  COMPARISON  FOR  TRUITT'S 
TRAJECTORY 

TIME  I  TRUITT  SKINTEMP  NACA 


(SECONDS)  (LB/FT^)  (LB/FT^)  (LB/FTM 
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m.  Tamparatura  Coapariaon  For  Truitt's  Trajectory 

Figure  6  and  Table  2  compare  Truitt's  freestream 
temperature  to  SKINTEMP's  freestream  and  reference  static 
temperature.  The  freestream  static  temperature  calculations 
compare  favorably.  Again,  the  divergence  in  values  can  be 
partially  attributed  to  differences  in  atmospheric  modelling. 

Local  temperatures  are  compared  with  the 
freestream  temperatures  to  illustrate  the  difference  of  using 
local  atmospheric  properties  instead  of  freestream  properties. 
SKINTEMP  uses  local  conditions  for  the  heat  transfer 
calculations  while  Truitt  utilizes  the  freestream  conditions. 
Temperatures  near  the  skin  of  a  hypersonic  vehicle  will  be 
much  higher  than  the  surrounding  freestream  air  due  to 
boundary  layer  interactions  (i.e.  skin  friction) .  Thus, 
calculations  using  freestream  properties  are  expected  to  be 
more  approximate  than  those  using  the  local  conditions. 


22 


TABLE  2.  TEMPERATURE  COMPARISON  FOR  TRUITT'S  TRAJECTORY. 


TIME 

(SECONDS) 


TRUITT 

FREESTREAM 

(“RANKINE) 


SKINTEMP 

FREESTREAM 

(“RANKINE) 


SKINTEMP 

LOCAL 

(°RANKINE) 


5492 


6291 
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b.  Reynolds  Number  Comparison  For  Truitt's 

Trajectory 

The  difference  in  theoretical  approaches  is  very 
apparent  when  comparing  the  Reynolds  nximbers  in  Figure  7  and 
Table  3.  Truitt's  Reynolds  number  calculations  are  very  close 
to  SKINTEMP's  initial  calculations,  but  again,  the  divergence 
in  freestream  atmospheric  properties  clearly  has  an  impact  on 
the  freestream  Reynolds  number.  The  only  significant 
consideration  is  that  this  difference  becomes  larger  as  the 
program  builds  on  data  from  the  previous  time  step. 

Figure  7  clearly  illustrates  that  the  reference 
state  Reynolds  number  is  significantly  smaller  than  the 
freestream  value.  Recall  that  the  reference  state  utilizes 
the  local  flow  field  values  at  the  edge  of  the  boundary  layer 
with  the  transport  properties  evaluated  at  the  reference 
temperature,  and  that  the  Reynolds  number  is  given  as 


Nl 


RE 


p*VX 

f  ' 


where 


(18) 


p  =  P/RT 
V  = 

and 

2. 27 £-8 

+  198.7 


(34) 

(35) 

(15) 
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The  Reynolds  number  is  Inversely  proportional  to  the 
temperature.  Use  of  the  local  values  for  temperature,  which 
are  significantly  larger  than  the  freestrecun  values  of 
temperature,  will  cause  the  Reynolds  number  to  decrease. 

Another  parameter  that  will  contribute  to  a 
difference  in  Reynolds  number  calculations  is  the  relationship 
used  to  determine  viscosity,  equation  15.  Truitt  did  not 
utilize  the  same  relationship  and  even  though  the  values  are 
of  the  same  magnitude,  the  difference  only  adds  to  the 
deviation  in  results. 

FREESTREAM  AND  REFERENCE  REYNOLDS  NUMBER  COMPARISON 


Figure  7.  Reynolds  Number  Comparison 
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TABLE  3.  REYNOLDS  NUMBER  COMPARISON  FOR  TRUITT'S  TRAJECTORY 


TIME 

(SECONDS) 

TRUITT  RE 
(  X  10’) 

SKINTEMP  RE 
FREESTREAM 
(X  10 

SKINTEMP  RE 
REFERENCE 
{  X  10’) 

1 

5.290 

5.039 

4.993 

2 

10.092 

5.051 

4.972 

3 

14.790 

14.081 

12.528 

4 

19.270 

19.329 

14.424 

5 

23.186 

22.908 

14.204 

6 

21.274 

22.693 

12.772 

7 

19.408 

21.189 

11.994 

8 

17.622 

19.188 

11.192 

9 

16.136 

15.839 

9.513 

10 

14.614 

15.066 

9.514 

11 

21.130 

18.134 

10.223 

12 

26.744 

20.591 

8.881 

13 

31.116 

21.301 

6.626 

14 

34.338 

20.262 

4.515 

15 

35.660 

19.013 

3.144 

16 

31.030 

17.679 

2.515 

17 

26.822 

13.764 

1.865 

18 

22.762 

12.733 

1.777 

19 

18.212 

10.984 

1.584 

20 

14.644 

8.447 

1.279 

21 

15.436 

6.840 

1.004 

22 

14.336 

6.380 

.770 

23 

12.356 

4.406 

.437 

24 

11.040 

3.513 

.298 

25 

8.280 

3.124 

.230 

26 

5.764 

1.724 

.117 

27 

4.212 

1.173 

.072 

28 

2.754 

1.171 

.071 

29 

1.902 

.906 

.049 

30 

.942 

.320 

.017 
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Skin  Temperature  (Rankine) 


c.  Skin  T&apejrature  Cotapariaon 

Figure  8  and  Table  4  ccMnpare  the  final  skin 
temperature  calculations  between  Truitt  and  SKINTBMP. 

TRUITT  SKIN  TEMPERATURE  PREDICTION  COMPARISON 


Time  (seconds) 


Figure  8.  Skin  Temperature  Comparison 
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TABLE  4.  SKIN  TEMPERATURE  COMPARISON  FOR  TRUITT'S  TRAJECTORY 


TIME 

TRUITT 

SKINTEMP 

(SECONDS) 

(“RANKINE) 

(“RANKINE) 

1 

520 

519 

2 

540 

522 

3 

570 

560 

4 

600 

633 

5 

660 

717 

6 

720 

772 

7 

700 

767 

8 

680 

742 

9 

660 

704 

10 

655 

677 

11 

665 

707 

12 

820 

850 

13 

920 

1091 

14 

1330 

1408 

15 

1460 

1761 

16 

1550 

1917 

17 

1620 

1845 

18 

1610 

1781 

19 

1550 

1711 

20 

1510 

1626 

21 

1560 

1704 

22 

1680 

1968 

23 

1740 

2132 

24 

2060 

2259 

25 

2150 

2375 

26 

2175 

2319 

27 

2250 

2297 

28 

2250 

2300 

29 

2250 

2310 

30 

2250 

2071 
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1.  Results  Prom  Truitt's  Trajectory 


The  graphs  and  tables  indicate  comparable  results  in 
the  two  different  methods  of  calculating  skin  temperature. 
There  are  numerous  reasons  for  the  slight  differences  in  the 
results.  The  discrepancy  in  standard  atmosphere  values  leads 
to  significant  differences  in  the  aerodynamic  parameters, 
which  cause  marked  deviations  in  the  heat  transfer 
coefficients,  and  the  adiabatic  wall  and  skin  temperatures. 

Another  major  difference  is  the  use  of  the  reference 
state  in  SKINTEMP  and  the  freestream  state  in  Truitt's  test 
case.  Using  the  reference  condition  instead  of  the  freestream 
state  is  considered  to  be  the  more  correct  model  of  what  is 
actually  occurring  aerodynamically.  One  aspect  of  this 
difference  is  manifested  in  the  Reynolds  number  calculation. 
Despite  the  significant  difference  in  Reynolds  numbers  used  to 
compute  the  Nusselt  number,  the  results  are  nearly  the  same. 
Recall  that  to  compute  the  Nusselt  number,  the  Reynolds  number 
is  raised  to  the  0.8  power.  Large  differences  in  Reynolds 
number  between  the  two  methods  are  minimized  in  this 
calculation,  allowing  the  skin  temperatures  to  more  closely 
correspond  with  one  another. 

Truitt  also  uses  varying  skin  material  specific  heats 
while  SKINTEMP  assumes  an  average  constant  specific  heat. 
This  does  not  cause  a  noticeable  change  in  skin  temperature 
calculations  but  does  account  for  a  few  degrees  of  difference. 
The  choice  of  time  constants  is  also  critical  for  stability 
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purposes.  Truitt  uses  a  time  step  of  two  seconds  while  the 
average  time  constant  step  generated  by  SKINTEMP  is  0.3 
seconds.  As  indicated  by  Truitt,  the  smaller  the  time 
constant,  the  more  accurate  the  results.  Truitt  used  the 
larger  time  constant  for  ease  in  calculations  because  they 
were  all  hand  computed.  With  all  other  variables  being  equal, 
the  calculation  of  skin  temperatures  with  the  more  simplified 
finite  difference  method  versus  the  Runge-Kutta  method  is 
expected  to  give  the  same  degree  of  accuracy,  provided  the 
time  step  is  small. 

While  the  Runge-Kutta  integration  method  is  an 
accurate  numerical  method,  Truitt  did  not  have  the  same 
computer  capabilities  that  are  in  use  today  and  the  current 
standard  atmosphere  tables  are  considered  to  be  more  accurate 
than  those  Truitt  used.  Even  though  the  SKINTEMP  results  are 
on  average  nine  percent  greater  than  Truittt's  results,  the 
skin  temperature  curves  are  similar  with  SKINTEMP  producing 
the  more  conservative  worst  case  results. 

A.  NASA  TEST  CASE 

NASA  AMES  provided  an  additional  test  case  for  SKINTEMP 
verification.  The  trajectory  is  typical  of  a  Single-Stage  - 
To-Orbit  (SSTO)  rocket  powered  vehicle  and  the  point  of 
interest  is  75  feet  aft  of  the  nose.  The  thickness  of  this 
point  of  interest  is  0.035  inches  and  the  material  at  this 
point  is  Tailorable  Advanced  Blanket  Insulation  (TABI) . 
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There  Is  a  body  angle  of  incidence  of  13  degrees  in  addition 
to  the  varying  angle  of  attack  experienced  through  flight. 
Figure  9  graphically  depicts  the  results  and  is  followed  by 
the  detailed  listing  of  skin  temperatures. 


NASA  SKIN  TEMPERATURE  PREDICTION  COMPARISON 


Figure  9 .  NASA  Skin  Ten^erature  Comparison 
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TABLE  5.  SKIN  TEMPERATURE  COMPARISON  FOR  NASA  AMES  TEST  CASE 


TIME 

(SECONDS) 


NASA 

("RANKINE) 


SKINTEMP 

(“RANKINE) 


6.82 


7.70 


8.71 


9.81 


10.76 


11.79 


12.85 


13.72 


14.72 


15.75 


16.75 


17.75 


18.75 


19.70 


20.78 


21.38 


2463 


2637 


2799 


2912 


3011 


3127 


3200 


3265 


3308 


3267 


3008 


3231 


2424 


2550 


2645 


2733 


2806 


2864 


2911 


2974 


3019 


3084 


2606 
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1.  Results  From  NASA  Trajectory 

Again,  the  results  from  the  two  skin  temperature 
calculation  methods  compare  favorably.  The  NASA  test  case  was 
also  based  on  freestream  conditions  with  surface  pressures 
given  by  tangent  wedge/cone  theory.  While  the  use  of 
freestream  conditions  is  not  theoretically  correct,  there  has 
been  significant  experimental  data  generated  to  justify  this 
assumption  [Ref.  14].  However,  these  assumptions  are  best 
used  for  upper  surfaces.  SKINTEMP  is  not  limited  to  this 
assumption  and  can  correctly  calculate  skin  temperatures  for 
upper  and  lower  surfaces. 
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V.  CGNCLUSIOHS  AMD  RBCOMMENDATIONS 

A.  CGNCLUSIONS 

The  SKINTEMP  con^uter  program  was  developed  as  a  tool  for 
use  in  the  conceptual  design  of  high  speed  missiles  and 
aircraft.  It  can  be  used  over  a  wide  range  of  trajectories, 
profiles  and  configurations  and  rapidly  yields  results  that 
enaible  the  user  to  easily  investigate  changes  in  material 
characteristics  and  flight  parameters  without  having  to  employ 
more  complex  and  cumbersome  calculations.  As  with  any 
computer  program,  the  results  can  only  be  as  good  as  the  data 
entered.  This  is  especially  true  of  SKINTEMP  because  valid 
results  can  only  be  expected  if  the  trajectory  is  accurately 
modeled. 

Validation  of  SKINTEMP  was  achieved  using  a  textbook  skin 
temperature  calculation  example  and  a  more  current  NASA  test 
case.  SKINTEMP  compared  with  the  NASA  and  Truitt  test  cases 
suggest  that  the  SKINTEMP  results  are  reasonable.  While  each 
of  the  methods  has  its  own  merits,  the  fundamental  theoretical 
base,  direct  approach  and  simplicity  of  SKINTEMP  make  it  ideal 
as  an  aerodynamic  and  astrodynamic  design  tool. 


B.  RECOMMENDATIONS 


SKINTEMP  currently  calculates  the  insulated  skin 
temperature  and  heat  flow  rate  at  only  one  surface  location. 
The  computer  code  could  be  expanded  to  perform  a  thermal  nodal 
analysis  that  could  predict  skin  teirperatures  along  an 
extended  surface.  Moreover,  ten^eratures  of  materials  within 
the  air/spacecraft  affected  by  the  skin  temperature  could  be 
calculated.  This  would  be  very  useful  in  the  design  process 
to  initially  identify  hot  spots  which,  especially  in 
hypersonic  flight,  could  be  crucial  to  the  overall  design. 
SKINTEMP  was  written  strictly  for  the  turbulent  flow  case 
because  this  generally  produces  worst  case  skin  temperatures. 
The  code  could  easily  be  modified  to  model  laminar  flow  if 
turbulent  flow  was  not  the  required  assumption. 

Because  the  skin  temperature  problem  is  quite  detailed  and 
the  inputs  are  extensive,  complete  trajectory  and 
corresponding  surface  skin  temperatures  profiles  are  difficult 
to  find  in  the  literature.  Additional  test  cases  could  be 
developed  for  SKINTEMP  for  further  validation  as  well  as  any 
experimental  data  that  is  available. 

The  important  objective  to  remember  in  any  significant 
code  modification  of  SKINTEMP  is  that  the  program  needs  to  run 
quickly  and  allow  the  user  to  modify  any  design  parameters. 
If  the  program  becomes  too  con^lex  and  time  consuming,  then 
the  usefulness  of  SKINTEMP  is  degraded. 
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APPENDIX  A.  SKZNTBMP  COMPUTER  CODE 


LT  Maty  L.  Cuulnga 
591 

This  progrsa  cslculstss  tcsnsisnt  alstloii  sarodynsalc  hasting 
psrsaatacs  (Adisbstlc  wall  taapacstutas.  local  hast  tcansCar 
coaffeiants,  and  Insulatad  akin  taaparaiuraa)  using  various 
basic  flight  paramatacs  as  Input. 


Varlabla  daclaratlons 

BODY  -  datarainas  shapa  of  body,  althar  wing  or 
axisyaaatrlc 

SNBEP  •  angla  of  wing  swaap  (dagcaas) 

TEMP  •  launch  vahlcla  taaparatuca  (Pahcanhalt) 

NUMOAT  ••  nuabar  of  lines  of  data 

TABLE  >  input  data  consisting  of  tlaa  (sac),  Mach  nuabar, 
of  attack  (dagraes),  and  halght  (feet) 

X  <■  distance  froa  leading  adga  to  point  of  Intarast  (feet) 

TUIN  •  thickness  of  Interest  point  (Inches) 

RAD  -  nose  radius  of  body  (ft) 

HATOENS  -  aatarial  density  (lba/ft'3) 

HATEniS  -  aatjarlal  aalsslvity 

HATSPHT  -  aatarial  speciflo  heat  (BTU/lba-r) 

TEHPRAO  -  radiation  taaparatura  (Pahrennheit) 

ANGLE  >  angla  froa  the  unit  surface  noraal  to  the 
freestraaa  (degrees) 

PRINTOUT  -  user  defined  print  tiae  (seconds) 
g  •  accaleration/welght  constant  ( Iba-f t/lbf-sac*2 ) 

LI  ••  taaparatura  lapse  rate  froa  0  -  36069  ft  (R/£t) 

L2  >  taaparatura  lapse  rata  froa  82021  -  1S4200  ft  (R/ft)  ' 
R  -  gas  constant  ( f t-lbf/slug-R) 

TH  -  thickness  of  intarast  point  (feet) 

Oansin  -  initial  density  at  sea  level  (slug/  ft‘3) 

Pin  -  initial  pressure  at  sea  level  (lbf/£t*2) 

Teapin  •  launch  vehicle  taaparatura  (Ranklne) 

T/ataT  -  frasstreaa  taaparatura  (Ranklna) 

P/atapr  ••  fraastraaa  prassura  (lbf/ft’2) 
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OENS/atad«nc  •  fr«tctcaaa  danalty  (slu^/  ft‘3) 

PrtisT  praicur*  at  tha  tropoaphaca  calling  (lb£/Ct*2) 
DansT  •  danalty  at  tha  tropoaphaca  calling  (alug/  Ct‘3) 
Ptl  -  local  stagnation  praasura  (lbC/ft*2) 

Pt  •  Craastraaa  stagnation  prassuta  (lbf/ft‘2) 

Tt  ••  Craastraaa  stagnation  taapatatuia  (Ranklna) 

Ma  ■  calculatad  fraaatraaa  Mach  nunbar 
Ml  •  Hach  nuabat  behind  shock  wava 
Hi  •  local  Mach  nuabar 
a  •  Iraaatraaa  sonic  velocity  (ft/sac) 

V  -  fcaastraaa  velocity  (ft/sac) 
q  ■  dynaalc  prassuta  ( slug/f t-sac' 2 ) 

gaaaa  -  ratio  of  specific  heat  at  constant  pressure  to 
specific  heat  at  constant  voluaa 

slgaa  •  Stafan-Boltsaann  constant  (BTU/hr-f t'2-R‘4 ) 

theta  ••  angle  of  attack  added  to  tha  nocaal  freastraaa 
angle  (radians) 

CpHax  -  Haxlaua  pressure  coefficient 
Pi  -  local  pressure  (lbf/ft‘2) 

Ttl  ••  local  stagnation  taaparatura  (Ranklna) 

Tl  >  local  taaparatura  (Ranklna) 

Tn  -  skin  teapertaura  (Ranklna) 

TAH  >  adiabatic  wall  taaparatura  (Ranklna) 
re  «  recovery  factor 
Pr  >  Prandtl  nuabar 
u  «•  viscosity  (slug/sac-ft) 
k  theraal  conductivity  (BTU/sec-f t-r) 

Cl  •  Sutharland’s  constant  ( slug/f t-sec-R' . 5 ) 

C2  •  Sutherland's  constant  (Ranklna) 

Cp  •  specific  heat  (BTU/lba-R) 

Tatar  -  rafaranca  taaparatura  (Ranklne) 
danstar  •  reference  density  (slug/ft*3) 
ustar  •  rafaranca  viscosity  (slug/sac-ft) 
al  -  local  sonic  velocity  (ft/sac) 
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VI  «  loc«l  velocity  (ft/««c) 

Ratter  -  rafaranca  Raynoldt  nusbac 

Cpttar  ••  rafaranca  tpaclflc  heat  (BTU/lbn-R) 

kttar  *  rafaranca  tharaal  conductivity  (BTU/aac-f t-F) 

Pretar  •  rafaranca  Prandtl  nuaber 

NUttar  •  rafaranca  Nuaeelt  number 

h  *  local  heatj  transfer  coefficient  (BTU/sec-f t‘2-F) 
w  >  specific  weight  (thf/ft*!) 

TBNPNEX  >  skin  teaperature  calculated  based  on  the 
previous  skin  temperature  (Rankine) 

Hr  -  time  differential  calculation  aid 

tlaa  tiae  step  (sec) 

suatiae  >  Interval  tiae  counter  (sec) 

runtiae  «  total  tiae  counter  (sec) 

printtlae  -  tiae  printing  counter  (sec) 

stagpt  -  stagnation  point  (ft) 

ANGLER  ••  Angle  of  attack  (radians) 

hsl  «  stagnation  enthalpy  (BTU/lba) 

hws  -  wall  enthalpy  (BTU/lba) 

hwS40  wall  enthalpy  at  540  R  (BTU/lba) 

qstag  *  stagnation  point  heating  rate  (BTU/sec) 

qheat  -  heat  transferred  (BTU/sec) 

dimension  TABLB( 100, 4 ) ,  T(IOO),  P(IOO),  DBNS(IOO) 
integer  BODY,  NUNDAT 

real  SWEEP,  TEMP,  TABLE,  X,  THIN,  RAD,  HATDENS,  HATENIS, 
C  HATSPHT,  TENPRAD,  ANGLE,  TH,  PRINTOUT,  g,  Ll ,  L2 ,  L3 ,  R, 
c  Densln,  Pin,  Tenpin,  T,  P,  DENS,  PressT,  DensT, 
c  atnpr,  atadens,  ataT,  Ptl,  Pt,  Tt,  HI,  a,  V,  q,  Ha, 
c  gaaaa,  sigaa,  theta,  pi,  printtiae, 

c  CpNax,  Pi,  Hi,  Ttl,  Tl,  Tn,  TAH,  re,  Pr,  u,  k.  Cl,  C2, 
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c  Cp.  Tatar,  danatar,  uatac,  al,  VI,  Raatac,  Cpatar, 
c  katar,  Pcatar,  NUatar,  ti,  w,  tbnpncx,  HR,  tiaa, 
c  atagpt,  ANGLBR,  hal,  hwa,  hwS40,  qatag,  qheat, 
c  tiaadlf,  auatiaa,  daltatiaa,  daltaa 


Input 


print  * 

print  *,  'If  tha  body  la  axisyaaetric,  enter  1.' 
print  *,  'If  the  body  ia  a  wing,  enter  2.' 

read  *,  BODY 

print  * 

if  (BODY  .EQ.  2)  then 

print  *,  'Hhat  la  the  wing  awaep?' 
read  *,  SHBBP 
print  * 
end  if 


print  *,  'Hhat  ia  tha  launch  teapecature  of  the  vehicle?  (r)' 
read  *,  TEHP 
print  * 

print  *,  'How  aany  data  pointa  will  be  entered  for  tiae,  Hach  nua 
c  her,  angle  of  attack,' 
print  'and  altitude?' 

read  *,  NUHOAT 

print  * 

print  *,  'Enter  tha  tiae  (a),  Hach  nuaber,  angle  of  attack  (degr 
c  eea),  and  altitude  (ft)  in  a  row/coluan  foraat.' 
print  *,  'Hit  return  after  every  entry.' 
print  * 

print  *,  ' ***HARNING***  Do  not  uae  any  data  linea  with  H  •  0' 

Do  20  ROH-1,  NUHDAT 
Do  10  COL-1,4 

read  *,  table(ROH,COL) 

10  continue 

20  continue 
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print  * 

print  *,  'What  ii  tha  dlitanca  from  the  tip  ot  the  vehicle  to  the 
c  point  of  intereit?  (ft)' 

read  *,  X 

print  * 

print  *,  'What  is  the  thickness  of  the  surface  at  the  point  of  In 
c  terest?  (in)' 

read  THIN 

print  * 

print  *,  'What  is  tha  nose  radius?  (ft)' 
read  *,  RAD 
print  * 

print  *,  'What  is  the  density  of  the  aiaterial?  (lbn/ft'3)' 
read  MATOENS 

print  * 

print  'What  is  the  ealssivity  of  the  material?' 
read  *,  HATEHIS 
print  *  I 

print  *,  'What  is  the  specific  n<.-.  j£  the  material  used?  (BTU/lb 
cm-r ) ' 

read  *,  HATSPHT 
print  * 

print  *,  'For  the  majority  of  the  flight  path,  if  the  surface  rad 
dates  towards  the  earth,  enter  60.' 

print  *,  'If  it  radiates  to  space,  enter  0.' 

read  «,  TENPRAO 

print  * 

print  *,  'What  is  the  angle  from  the  unit  surface  normal  to  the 
cfreestream?  (degrees)' 

read  *,  ANGLE 

print  * 

print  *,  'What  time  increment  would  you  like  the  skin  temperature 
cs  printed?  (seconds)' 

read  *,  PRINTOUT 

print  * 
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50 


gaMa  >  1.4 

algaa  -  .1714E-8/3600.0 
pi  -  3.1415927 
TH  -  THIN/12.0 
auatlaa  •  0. 
deltatiaa  -  0. 
pclnttiaa  -  printout 

if  (TEHPRAD  .EQ.  60)  then 

TEHPRAO  -  TEHPRAD  +  460. 

end  if 


do  100  n  -  1,  NUHOAT 

if  (n  .EQ.  1)  then 
Tn  -  Tenpin 
«a  -  TABLB(n,2) 
atnpc  >  P(n) 
atadens  -  Denain) 
ataT  -  T(n) 

if  (BODY  .EQ.  2)  then 

Ma  -  Hp  *  cos(sweep  *  pi/180.) 
end  if 

tiaedif  -  (TABl.E{n+l,l)  -  TABLB(n.l)) 
GOTO  60 
end  if 

deltaa  -  (TABLE(n-t-l,2)  -  TABLE(n,2)) 
tiaedif  -  (TABLB(n+l,l)  -  TABLE(n,l)) 

if  (auatine  .GT.  0.)  then 
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AtBo«ph«rlc  data  calculatlona 


a 

a 


g  -  32.174 
LI  -  -.0035665 
L2  -  .001645 
L3  -  -.002469 
R  -  1716. 

Oansln  -  .002378 
Pin  -  2116.0 
Teapln  -  TEMP  *  460.0 

do  30  HOT  «  1,  NUNOAT 

if  (TABLE(UGT,4)  .LB.  36089.0)  than 
Tcopoaphere 

T(HOT)  -  Taapln  *  Ll  *  TABLE(HGT,4) 

P(HGT)  -  Pin  *  (t(HGT)/Te«pln)**(-g/(Ll  *  R  )) 

OBNS(HGT)  -  Oensin  *  < (T(HGT)/TeBpln)**-(g/(Ll  *  R) 

c  ♦  1)) 

goto  30 

end  If 

Stratosphere 

If  (TABLB(HGT,4)  .GB.  82021)  then 

if  (TABLE(HGT,4)  .GE.  173885)  then 

T(HGT)  •  509  *  L3  * ( TABLE( HGT, 4 )  -  173085) 

P(HGT)  -  1.5  *  (T(HGT)/509)**(-g/(L3  *  R  )) 

OENS(HGT)  •  .0000014  *  (  (  T(  iIGT)/509  )  *  * 

C  -((g/(L3  *  R))+  1.0)) 

goto  30 

and  If 

if  (TABLE(RGT«4)  .GB.  154200)  then 
T(HGT)  -  509 

PressT  -  2.6419  / 
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OeniT  -  .00000303 

P(liGT)  -  Pr**ST  *  BXP  (-((g  /  («<  *  T(HGT)))  • 
C  (TABLB(HGT,4)  -  154200.))) 

DEN8(UGT)  -  OBnsT  *  P(HGT)  /  PCessT 

goto  30 

•nd  1C 

T(HGT)  -  389.9  t  L2  *  ( TABLE ( UGT , 4 )  -  82021.) 

P(HGT)  -  55.8  *  (T(HGT)/389.97)**(-g/(L2  *  R  )) 

OENS(HGT)  -  .0000789  *  ( ( T( UGT)/389 . 97 ) *  * 

C  -((g/(L2  •  R))+  l.Oi) 

goto  30 

and  iC 

T(HGT)  -  -67  <  460.0 

Pr«BsT  ■>  476 
DansT  -  .000710 

P(HGT)  -  PraaaT  *  EXP  (-((g  /  (R  *  T(HGT)))  *  ( 
c  TABLB(HGT«4)  -  36089.))) 

DENS(8GT)  -  OansT  *  PlHGT)  /  PraasT 
30  contlnua 


Haadar 


print  * 
print  * 
print  * 
print  * 

print  40,  'Tina',  'Adiabatic  Hall',  'Haat  TransCar', 
c'Skin  Teap',  'Haat  TransCar  Rata  ' 

print  45,  '(sac)',  '  Taaparatura  (R)',  ' CoeCCiclant ' ,  '  (R)', 

c  '  (BTU/sac-ona  sq  Ct)' 


print  * 


40 

roraat(2X, 

A4, 

2X, 

A16, 

IX, 

A16,  3X, 

A9, 

5X,  A21) 

45 

Por«at(2X, 

A5, 

2X, 

A16, 

4X, 

All,  6X, 

A4, 

8X,  A21) 

*  Subsonio  vs  supatsonic  calculations 
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n*  -  luktla*  *  d*lt«a/( tlBcdlf )  *  TABLED n. 2) 


c 


atMpc  ••  P(n)  4  sumtlaa  *  (P(n-fl)  -  P(n)  )/tlM«dl£ 

atadans  -  Dena(n)  4  suMtlae  *  (Dens(n4l)  -  Dan6(n)) 
/tlMdlf 

ataT  -  T(n)  4  suatlae*(T(n4l)-T(n) )/tiBedlf 

if  (BODY  .EQ.  2)  then 

Na  -  Ha  4  costsweep  *  pi/180.) 
end  if 


end  if 

if  (suatiae  .EQ.  0  )  then 

Ha  •  TABLE(n,2) 
atapr  •  P(n) 
atadens  -  Oena(n) 
ataT  •  T(n) 
if  (BODY  .EQ.  2)  then 

Ha-Ha  *  cos(6weep  *  pi/180.) 
end  if 

end  if 
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c 


c 


c 

c 

c 


if  (Ha  .LE.  1.0)  then 

Ptl  -  atapr  *  (1.0  4  ( ( gamaa-l . 0 )/2 . 0 )  *  na*42)** 
(gaaaa/(gaaBa-1.0) ) 

else 

Pt  -  atapr  •  (1.0  4  ( ( gaana-1 . 0 )/2 . 0 )  *  Ha**2)** 

( gaaaa/( gaaaa-1 . 0 ) ) 

Tt  -  ataT  *  (1.0  4  ( ( gaana-1 . 0 )/2 . 0 )  *  Ha* *2) 

Ptl  -  Pt  *  ( ( (gaaaa4l.0)/2  *  Ha**2)/(1  4  (gamma 
-1.0)/2  *  na**2 ))**( gamma/( gamma-1 . 0) )  * 
(2.0*gaaBa/(gamma4l.0)  *  Ha**2  -(gamma-1.0)/ 
( gaBma4l .0) )**(!/( 1-gamma ) ) 


end  if 

Adiabatic  wall  temperature  calculations 


4.6 


a  -  (gaaaa  *  R  *  atBT)**.5 
V  -  a  *  Ma 

q  -  .5  *  atadans  *  V**2 

theta  -  (ANGLE  *  TABLE(n,3))  *  pi/iao. 

CpMax  -  (Ptl  -  atapr)/q 

Pi  •>  q*CpHax  *  (co8(theta)  )**2  +  atapr 

Hi  •  ( ( ( Ptl/Pl ) **( (gaaaa-1 .0 )/gaaaa)-l . 0)  *  (2.0/ 
c  (gaaaa-1.0) ) )**.5 

Ttl  -  ataT  *  (1.0  +  ( gaaaa-1 . 0 )/2 . 0  *  Ma**2) 

T1  -  Ttl/(1.0  +  ( gaaaa-1. 0)/2.0  *  Ml**2) 

k  -  (.0141  *  (717.0/(T1  +  225.))  *  (Tl/492.0) 
c  **1.5)/3600.0 

Cp  -  .219  +  .342E-4  •  Tl  -  .293E-8  *  Tl**2 

Cl  -  2.27E-e 

C2  -  198.7 

U  -  Cl  *  Tl**1.5/(T1  +  C2) 

Pr  -  Cp  *  u  *  g  /k 
ce  ■  Pc**<l. 0/3.0) 

TAW  -  T1  *  (1.0-i^re  *  (gaBaa-1.0)/2.0  *  Hi* *2) 


Heat  transfer  coefficient  calculations 


Tstar  -  Tl  +  .5  *  (Tn-Tl)  +  .22  *  (TAW-Tl) 

denstar  >i  Pl/(R  *  Tstar) 

ustar  -  (Cl  *  Tstar»*1.5)/(Tstar  +  C2) 

al  -  (gaaaa  *  B  *  Tl)**.5 

VI  -  al  *  Hi 

Restar  ••  denstar  *  vl  *  x/ustar 

Cpstar  -  .219  +  .342E-4  *  Tstar  -  .293E-B  *  Tstar 
c  **2 

kstar  -  (.0141  *  ( 717 . 0/( Tstar  +  225.))  *  (Tstar/ 
c  492.0)**1.5)/3600.0 

Prstar  •  Cpstar  *  g  *  ustar/kstar 
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c 


c 

c 


if  (Rastar  .LK.  1.B7)  than 

NUatar  .0292  *  (Restar**. 6)  *  (Prstar 
**(1. 0/3.0)) 

else 

NUstar  •  .185  *  (Restar)  *  (Prstar)** 
(i.0/3.0)/( ( logl0(Re6tar ) )** 
2.584) 

end  i  £ 

h  -  NUstar  *  kstar/x 


Skin  temperature 


calculations 


W  -  MATDENS  /  9 

Hr  -  MATEMIS  *  Sigma  *  (Tn**3  -t-  Tn**2  *  TEMP 
C  RAO  *  TEHPRAD**2  *  Tn  TEHPRAO**3) 

time  -  .25  *  MATSPHT  *  w  *  g  *  th  /(hr  +  h) 

TEHPNEX  -  time*(((h  *  1  *(TAH-Tn)) 
c  -  sigma  *  MATEMIS  *  (Tn**4-Tr**4 ) )/ 

C  (MATSPHT  *w*g*TH*l))+Tn 


I 


Stagnation  point  heat  transfer  calculation 


ANGLER  -  ANGLE  *  pi/180. 

hsl  -  Cp  *  Ttl 

hws  -  Cp  *  Tn 

hw540  •  Cp  *  540. 

stagpt  ••  cos(70  *  (pl/100.)) 

if  (ANGLER  .LE.  Stagpt)  then 

qstag  ••  17600. 0/SQRT( RAD)  *  S0RT( atmdens/Densin )  *  (V/ 
c  26000. 0)**3. 15  *  ( (hEl-hws)/(hEl-hw540) ) 

if  (BODY  .EQ.  1)  then 

qheat  -  qstag  *  (cos(ANGLER) ** ( 3 .0/2 .0) ) 
else 

qheat  -  qstag  *  ( cos(ANGLER) ** ( 3 . 0/2 . C ) )/SQRT( 2 . 0 ) 
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•nd  iC 


•  Isa 


Uaat  transfac  calculation 


qhaat  •!>*!*  (TAH-Tn)  -  sigaa  * 

C  NATEHIS  *  (Tn**4-TEHPBAD**4) 

•nd  If 

Tn  •  Taapnax 

suatiaa  •  suatlaa  +  tia« 
runtiaa  >  runtlaa  *  tla« 

Tf  •  Taapnax 

If  (suatiaa  .GE.  tiaadif)  than 
daltatiaa  •  suatiaa  -  tiaadif 

if  (runtiaa  .GE.  printtiaa)  then 

Print  80,  printtiaa,  TAH,  h,  Tf,  qhaat 
HRITE(3,90)  printtiaa,  TEMPNEX 
printtiaa  «  (printtiaa  printout) 

•nd  if 
suatiaa  ••  0 

goto  100  ^ 

else 

if  (n  .GE.  nuadat)  then 
goto  100 
end  if 
GOTO  50 
end  if 


Output 


* 

* 

* 
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80 

rorMt(  OX,  r7.2,  SX, 

r9.2,  lox,  re. 6,  sx,  r9.2,  i2x.  rs.a) 

90 

roraatdX,  r7.2.  lOX, 

r9.3) 

100 

continue 

print  * 
print  * 
print  * 

end 


I 


50 


APPENDIX  B.  EXAMPLE  SKINTEMP  INPUT  AND  OUTPUT 


Wkiit  111  till*  Imifirli  t  i*ipfi«>riitiif  I*  of  the  ttuhlcle?  tri 

Now  iMiit  d-itt  wilt  h«  •intiitRil  let  t Ino,  HjirIi  nimi  but,  anqle  of  attacli. 

and  altitiifin? 

tntni  the  tlin*  (at «  t<ach  nimiliRt ,  nnifle  e{  attach  (iloir  ecal ,  and  altitude  ((tl 
In  a  tew/ceiuian  letmat. 

Nit  return  after  every  entry. 

*  • ‘MANMIMfS*  *  *  no  not  ii*e  any  data  line*  with  M  «  0 

What  |»  the  diatance  from  the  tip  of  the  vehicle  to  the  point  of  Interest?  (ftl 
What  Is  the  thickness  of  the  surface  at  the  point  of  In  tereat?  finl 
What  la  the  nose  ladlua?  <(tl  * 

What  Is  the  density  of  the  natertal?  (Ibni/lt*3} 

What  Is  the  emlssivlty  of  the  material? 

Whst  Is  the  specific  heat  of  the  material  used?  fOTO/lhm-rl 

loi  the  majority  of  the  lllyht  path.  If  the  surface  rndlatea  towards  the  earth, 
enter  *0. 

If  It  radiates  to  apace,  enter  0. 

What  Is  the  nn^le  from  the  unit  surface  normal  to  the  freestream?  Ideyreesl 
What  time  Increment  would  you  like  the  akin  temrreratures  printed?  fsecondsi 

SKINl'EMP  data  entry  prompts  ~ 
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1 

a 

i7 

1 

.3( 

0 

2^0 

2 

.US 

0 

900 

3 

l.Ot 

0 

390 

4 

1.449 

0 

1000 

9 

l.SK 

0 

1900 

( 

1.339 

0 

2000 

3 

l.«S 

0 

3000 

t 

1.930 

0 

9000 

t 

1.48 

0 

*900 

10 

1.3*3 

0 

10000 

11 

2.103 

0 

13000 

12 

2.83 

0 

19000 

13 

3.939 

0 

18000 

14 

4.348 

0 

22000 

19 

9.153 

0 

25000 

1C 

9.0*5 

0 

23000 

13 

5.03 

0 


entered  usine  the  file  etil 


litt  te.<!t  ca.se> 


» 


34000 

to 

4.tSt 

0 

37000 

to 

4.004 

0 

40000 

20 

4.(49 

0 

4(000 

21 

s.oo 

0 

soooo 

22 

(.951 

0 

S2000 

23 

0.103 

0 

(2000 

24 

9.294 

0 

((000 

29 

to. 373 
0 

(9000 

2( 

tt.92t 

0 

09000 

27 

i2.(7 

0 

09900 

20 

33.017 

0 

99000 

29 

14.43 

0 

117000 

30 
15 
0 

129000 

31 
IS 
0 

125000 

32 
IS 
0 

125000 

20 

.005004 

.05 

530 

.9 

.121 

0 

90 

1 


Data  entered  using  the  file  entry  method  (Truitt  test  case,  continued) 
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T  i  i»M ) 

A(.liAb»tlc  W.Tll 

IlcAt  Tranalet 

Skin  Temp 

Heat  Trartsfer  Rate 

(SPCl 

Tt^mfierAture  (P) 

Coefficient 

(R1 

(BTO/sec-one  sq  ftl 

1.00 

529.30 

0.008575 

522.02 

0.041 

2.00 

609.95 

0.019815 

560.91 

1.244 

3,00 

698.44 

0.024123 

633.11 

2.018 

4.00 

793.95 

0.025823 

717.28 

2.503 

5.00 

794.17 

0.024537 

772.03 

0.526 

r.  .00 

760.86 

0.023045 

767.06 

-0.392 

7.00 

725.14 

0.021186 

742.17 

-0.660 

:>  .00 

«74.31 

0.017698 

104.24 

-0.8S4 

•J.OO 

656.61 

0.017273 

671.12 

-0.611 

10.00 

792.45 

0.019156 

707.28 

2.048 

n  .00 

1042.91 

0.019434 

850.20 

4.754 

12.00 

1397.34 

0.017993 

1091.53 

6.712 

13.00 

1859.94 

0.0IS438 

1408.19 

7.561 

H.OO 

2451.81 

0.013316 

1761.73 

7.877 

15.00 

2457.28 

0.011619 

1917.23 

0.815 

16.00 

2257.59 

0.00(188 

1845.07 

-2.340 

17.00 

2202.84 

0.008247 

1781.19 

-1.444 

18.00 

2096,57 

0.007787 

1711.37 

-1.503 

13.00 

1979.69 

0.005898 

1626.52 

-1.607 

70.00 

2740.36 

0.005182 

1704.54 

3.274 

71.00 

3870.34 

0.004645 

1967.57 

5.554 

72.00 

5010.18 

0.003351 

2131.51 

2.475 

73.00 

6573.34 

0.002756 

2259.56 

2.906 

74.00 

8142.86 

0.002458 

2374.87 

2.701 

75.00 

9968.70 

0.001562 

2319.48 

-2.354 

2fi.OO 

12715.81 

0  001156 

2297.95 

0.757 

O 

o 

12728.18 

0.001152 

2300.75 

0.071 

28.00 

15484.49 

0.000928 

2310.05 

0.229 

18435.02 

0.000420 

2071,97 

-5.430 

30.00 

70525.59 

0.000321 

1952.78 

-1.970 

11  .00 

70525.59 

0.000323 

1936,63 

-0.227 

ir.oo 

70525.59 

1 

0.000324 

1935,69 

-0.013 

SKINTEMP  ou4>ut  (Truitt  test  case) 
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